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£/ Abstract 

pq , In this article we have calculated the spin-independent cross section of nucleon- 

dark matter scattering process at loop level, which is relevant to dark matter direct 
detection. Paying particular attention to the scattering of gluon with dark mat- 
ter, which contributes as leading order in the perturbation, we have systematically 
evaluated loop diagrams with tracking the characteristic loop momentum which 

C^~ ■ dominates in the loops. Here loop diagrams whose typical loop momentum scales 

are the masses of quarks and other heavier particles are separately presented. Then, 
we have properly taken into account each contribution to give the cross section. We 
assume that the dark matter is pure bino or wino in the supersymmetric models. 
The application to other models is straightforward. 



1 Introduction 

The existence of nonbaryonic dark matter has been established by cosmological obser- 
vations [I]. Weakly- interacting massive particles (WIMPs) are attractive candidates for 
dark matter, and many models are proposed to predict WIMPs. The lightest neutralino 
in the minimal supersymmetric standard model (MSSM) has been extensively studied 
among them. 

The direct detection experiments of WIMP dark matter are one of the methods for 
probing the nature of the dark matter. In the experiments, one searches for the signatures 
of WIMP-nucleon scattering. Many experiments have searched for the dark matter signals, 
and their sensitivities have been improved. The CDMS II experiment reported the final 
result for the five-tower WIMP search last year [2J. They observed two events in the 
dark matter signal region for the detector; their data set the upper limit on the spin- 
independent WIMP-nucleon elastic scattering cross section (donated as <7g[) of 3.8 x 
lO" 44 cm 2 for a WIMP of mass 70 GeV. This year the XENON100 experiment also 
started, and has already reported the first result [3]; they have given the upper limit as 
a^ < 3.4 x 10" 44 cm 2 for WIMPs of mass 55 GeV. In addition, the XMASS experiment will 
start soon [I]. Prospects for sensitivities of these new experiments are around 10 -45 cm 2 , 
and they cover broad parameter space of the neutralino dark matter scenario in the 
MSSM. 

In order to study the nature of dark matter based on those latest experiments, we 
need to evaluate the WIMP-nucleon elastic scattering cross section precisely. The WIMP- 
nucleon elastic scattering is induced by effective interactions of the WIMP with quarks 
or gluon at parton level. The effective interaction with the gluon is typically of higher 
order of the QCD coupling constant than those with quarks in many models, since the 
WIMP interacts with the gluon via loop diagrams including quarks and/or other colored 
particles. However, the effective scalar interaction with the gluon may make the leading 
contribution in the spin-independent (SI) cross section. This comes from the fact that 
the nucleon mass is dominated by gluon contribution. 

The effective interaction of the WIMP with the gluon is generated by loop diagrams. 
The loop momentum which dominates the integration is characterized by the mass scale 
of either heavy particles (such as WIMP mass) or quarks. We call the former (the latter) 
the "short-distance" ("long-distance") contribution. The long-distance contribution to 
the effective scalar operator of the WIMP with gluon is approximately evaluated from the 
effective scalar operators with quarks by using the trace anomaly of QCD. On the other 
hand, the short- distance contribution needs to be explicitly calculated. Drees and Nojiri 
calculated the one-loop box diagrams for gluon-neutralino scattering in the MSSM, and 
included the short- distance contribution to the SI cross section [5] . Although the SI cross 
section is evaluated in various models, the short-distance contribution is not included in 
many papers. 

In this paper we show a method to include the short- and long-distance contributions 
to the SI cross section systematically. We use the Fock-Schwinger gauge for the gluon 
field. This is frequently used in the QCD sum rules in order to evaluate the operator 



product expansions [B]. It is found that this gauge makes the calculation for the gluon 
contribution to the SI cross section much more transparent. In addition, we separately 
calculate short- and long-distance contributions, which makes it possible to add each 
contribution properly to the cross section. 

So as to study the details of gluon-neutralino scattering, we consider two cases where 
the lightest neutralino is pure bino and wino in the framework of the MSSM. For the 
bino-like neutralino case, the neutralino-quark scattering is induced by squark exchange at 
tree-level, and the neutralino-gluon scattering is generated at one-loop level. Calculating 
all the relevant diagrams at leading order, consequently we have found that the long- 
distance contribution of light quarks is partially regarded as the short-distance one in 
the calculation as Ref. [3]. This correction changes the SI cross section slightly in the 
MSSM. On the other hand, for the wino-like neutralino case, when superparticles except 
for the SU(2)l gauginos and the heavier Higgs bosons are decoupled in the MSSM, the 
neutralino-nucleon scattering process is dominated by weak gauge boson loop diagrams 
[7] . We evaluated the effective interactions of the neutralino with quarks at one-loop level 
and those with gluon at two-loop level. Here we use the Fock-Schwinger gauge for the 
gluon field, again; the result has been published in Ref. [S] • In this paper, we show detail 
of the evaluation. 

Our method to evaluate the effective interaction of the WIMP with the gluon is ap- 
plicable to other models of dark matter, such as universal extra-dimension scenario. Such 
analysis will be given elsewhere [9]. 

This article is organized as follows. In Sec. [2] we summarize the effective couplings of 
the WIMP with quarks and the gluon. We assume that the WIMP is Majorana fermion. 
Here we explain how to evaluate the effective interaction of the WIMP with the gluon. In 
Sec. [3j assuming that the neutralino is bino-like, we evaluate the effective couplings of the 
neutralino with quarks and the gluon by using the Fock-Schwinger gauge for the gluon 
field. Then we consider the case that the neutralino is wino-like where the neutralino- 
nucleon scattering is generated at loop-diagrams due to the weak interaction in Sec. HI 
Finally Sec. [5J is devoted to the conclusion. In the Appendix, we show the propagators 
of scalar and fermion in a gluon background (Fock-Schwinger gauge) and other useful 
formulas for evaluation of effective interaction of the WIMP with the gluon. 

2 Effective interaction for WIMP-nucleon scattering 

First, we summarize the effective interactions of the WIMP with light quarks (q = u, d, s) 
and gluon, which are relevant to the WIMP-nucleon spin-independent scattering. They 
are given as follows, 

£cff = £ £ eff + £ efT (1) 

q=u,d,s 



where 

(1) (2) 

Cf = fGXXG^G^ + ^^id^YxO^ + ^^id^idnxO^. (3) 

X is the WIMP which we assume to be Majorana fermion as it is mentioned in the 
Introduction. M and m q are masses for the WIMP and quarks, respectively. The second 
and third terms in C e q s and £^ ff depend on the twist-2 operators (traceless parts of the 
energy-momentum tensor) for quarks and gluon, respectively, 
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°lv = o?M D iH» + A,7*« - o9w P)Q > 



0% = (^G^ + \ 9lv Gty?*\ . (4) 

Here the covariant derivative is defined as D^ = d^ + ig s A^T a (g s , T a and A® are the 
SU(3)c coupling constant and generator, and the gluon field, respectively), and G" is 
the field strength tensor of gluon. 

The SI cross section of the WIMP with target nuclei T is expressed compactly in terms 
of the SI coupling of the neutralino with nucleon f N (N = p, n) [10]; 
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*si = -^r^J \n P f P + n nfn \ , (5) 

where my is the mass of target nucleus, and n p and n n are proton and neutron numbers 
in the target nucleus, respectively. The SI coupling of the neutralino with nucleon is given 
by the coefficients and matrix elements of the effective operators in C q s and C e g s ad_| 

f N /m N = ^ f q f Tq + J2 ^(2) + g(2))(s«+4 2 )) 

q=u,d,s q=u,d,s,c,b 

~ £/"/« + |g< 2 ) (sS'+sg 1 ) ■ (6) 

The matrix elements of the effective operators are expressed by using nucleon mass as 

(N\m q qq\N)/m N = f Tq , 

1 - 2_^ fcq = JtG , 

u,d,s 

(N(p)\Ol\N(p)) = J-(p^-i m ^)(g(2) + 5(2)), 

777 ^y 4 

(N(p)\0^\N(p)) = J-( p ^-i m ^)G(2). (7) 

777 ^y 4 



1 We have included contributions from effective interactions with twist-2 operators for charm and 
bottom quarks, g£ and g b % (i = 1,2), in Eq. ©. See later discussion. 



In the matrix elements of twist-2 operators, q(2), g(2) and G(2) are the second moments 
of the parton distribution functions (PDFs) of quark, antiquark and gluon, respectively, 

q{2) + q(2) = / dx x [q(x) + q{x)\ , 
Jo 

G{2) = dx x g{x) . (8) 

Jo 

For the numerical calculation, we describe the input parameters for the hadronic ma- 
trix elements that we use in this article. The mass fractions of light quarks, /y ? , are 
expressed as [UJ 

fru = -^(l±fl^, 

m u + m d m N 

hi = ; (i=fO — , 

m u + m d m N 

hs = ; y — (9) 



for proton/neutron. Here, 



m u + m d - 

v-kN = {p\uu + dd\p) 



2(p\ss\p) 
(p\uu + dd\p) 
= {p\uu - dd\p) 
(p\uu + dd\p) 

The recent lattice simulation predicts that cr^jv = (53 ± 2(stat)l7 1 (syst)) MeV and y = 
0.030 ±0.016(stat)+°;°^(extrap)+^(m s ) pj. Combined with f = 0.132 ±0.035 PUCE], 
we get f Tu = 0.023, f Td = 0.032, and / Ts =0.020 for the proton, and / Tu =0.017, f Td = 0.041, 
and hs— 0.020 for the neutron. In our numerical evaluation, we use the center value for 
each parameters. 

The second moments of the PDFs of quark, antiquark and gluon are scale-dependent, 
and are mixed with each other once the QCD radiative corrections are included. Thus, the 
coefficients of the terms with twist-2 operators become scale-independent after they are 
multiplied by the second moments of the PDFs [5] . We use the second moments for the 
PDFs at the scale of Z boson mass, and include bottom and charm quark contributions, as 
given in Eq. fl6]). The second moments are evaluated from the CTEQ parton distribution 
as 

G{2) = 0.48, 

«(2) = 0.22, 0(2) = 0.034 , 

d(2) = 0.11, J(2) = 0.036, 

s(2) = 0.026, s(2) = 0.026 , 

c(2) = 0.019, c(2) = 0.019 , 

6(2) = 0.012, 6(2) = 0.012, (11) 



for the proton. Those for the neutron are given by an exchange of up and down quarks. 

For the calculation of f^ at leading order, we give some remarks. The effective inter- 
actions of the WIMP with the gluon are of higher order of a s (= g^/Air) than those of the 
light quarks in many models including the MSSM. However, it is found in Eq. (Q that the 
effective scalar coupling of the WIMP with the gluon fa gives the leading contribution to 
the cross section even if it is suppressed by one-loop factor compared with those of light 
quarks. On the other hand, the contributions with the twist-2 operators of the gluon are 
of the next-leading as g G and g G are suppressed by a s . Then, we ignore them in this 
paper. 

Finally we explain the evaluation of fa- In general, the effective coupling of \-g 
scattering is induced by loop diagrams in which virtual quarks and other heavy particles 
run (like squarks in the MSSM). Consequently, the integral momentum around not only 
the quark mass scale but also the heavy particle one contributes. (As we described in the 
Introduction, we call them as "long-" and "short-distance" contributions, respectively.) 
Thus, in order to take into account both contributions accurately for general cases, we 
must calculate loop diagrams explicitly. So as to calculate diagrams including the quark 
systematically, we separate the two contributions as 

IgU = fdq + fdq ■ (12) 

In the long-distance contribution, there is another approximated way to evaluate the 
heavy quark (Q = c, b, t) contribution, which is valid in the limit of large mass of heavy 
particles. After integrating out heavier particles than heavy quarks, we write down the 
effective Lagrangian for heavy quarks: Cq = JqvtiqxxQQ + •••□• We calculate one-loop 
triangle diagrams in which heavy quarks rotate to emit two gluons, and geto 

/gIq = ~^Z^Q ■ ( 13 ) 



(See Appendix for detail.) Thus, we expect that 



a s 

12V 



fo\^ - ~fq , (14) 



in the limit of the large mass of heavy particles. (This can be utilized to check explicit 
loop calculation. See the following sections.) In the explicit loop calculation, we observe 
that the momentum with the scale of quark mass dominates the integral to give a factor 
1/rriQ (which is canceled in JqviiqxxQQ\ 



2 We used the equation of motion for the quark to reduce the effective Lagrangian. Though the 
heavy quarks are in loops in the evaluation of the effective interaction of the WIMP with the gluon, 
it is justified to use it. See Ref. [13] and references in it. We also explicitly checked that the effective 
interaction XX Qi PQ leads to the same result for fa as tuqxx QQ- 

3 This result is also obtained in another way. It is pointed out in Ref. [SJ that the effective interactions 
which come from heavy quark loops are approximately calculated in the effective interaction of the WIMP 
with heavy quarks, using the trace anomaly of energy-momentum tensor in QCD 16 . When the virtual 
heavy quark is the only particle in the loop diagram, the integrated-out heavy quark converts to the 
gluon operator as rriQQQ — > — ^G" G '"', which gives the same result as Eq. (TU 



For light quark contributions in long-distance loops, on the contrary, we must not 
include them in the same manner as heavy quark ones when we evaluate fa- As we saw 
in the above calculation, the loop momentum is dominated by quark mass; thus, the 
propagators for the soft quark whose momentum and mass are smaller than the QCD 
scale are determined entirely by confinement dynamics. The corresponding effect should 
be included in the matrix element (N\qq\N) jH]. That is why the light quark loops must 
not be counted in the evaluation of fa- 

Being aware of this fact, we sum up for quarks properly and get 

/« = £/!%+ E $V as) 

q=a\\ Q=c,b,t 

In addition, it is pointed out [17] that the long-distance contributions to fa suffer from 
large QCD correction. The corrections to Eq. (TT3]) are given as 

OL 

Ig\q = -tttCqIq > (16) 

LllX 

where cq = 1 + lla s {mQ)/<iii. We take c c = 1.32 and q, = 1.19 for a s {m z ) = 0.118 {m z 
is Z boson mass), and q = 1 in our numerical calculation. Then, rewriting Eq. ( TT51) with 
the QCD correction, we finally get 

5=all Q=c,b,t 

In the case where the WIMP interacts with the Higgs boson, the short- distance con- 
tribution to the effective interaction of the WIMP with the gluon is not generated at the 
leading order. Then, we can easily evaluate the SI cross section without explicit calculation 
of loop diagrams by using Eq. (fl6l) . However, it is not the case for the bino-like neutralino 
dark matter. In the next section, we consider the case where the lightest neutralino is 
almost pure bino in the MSSM. The bino interacts with quarks and squarks. We will see 
that the squark-quark loop diagrams give both short- and long-distance contributions. 

3 Bino-like neutralino dark matter 

The lightest neutralino interacts with quarks and gluon in diagrams with Higgs boson 
exchange or squark exchange. When the Higgsino-gaugino mixing is not negligible, the 
Higgs boson exchange tends to dominate in the SI x~N scattering. On the other hand, 
when the mixing is negligible and squark masses are comparable to the neutralino mass, 
the squark exchange dominates over the Higgs boson exchange. In this section, we take a 
limit \x ^> M {p, is the Higgsino mass parameter in superpotential) and consider the case 
where the lightest neutralino is almost pure bino. The extension to the general cases is 
straightforward. 

Before going to the numerical calculation, we show the couplings of the effective La- 
grangian which are given by the explicit calculation of tree- and loop-level scattering 





Figure 1: Tree-level diagrams of squark (q) exchange to generate interaction of WIMP 
and light quark. 



amplitudes. The following calculation is applicable to the models in which the WIMP 
interacts with quarks and new colored scalars. (For the time being, we do not distinguish 
a light or heavy quark for "g" unless it is expressly described, and provide the analytic 
results of the effective couplings.) 

First, we simply parametrize the neutralino interaction with quark and squark q as 

C = q(a q + b q >y 5 )xq + h.c. (18) 

Integrating out squarks at tree level (the corresponding diagram is shown in Fig. [T]), the 
effective interaction with light quarks in Eq. ([2]) is generated as 



/„ = 


M a 2 q + b 2 q 
(mf - M 2 ) 2 8 


(m§- 


1 a 2 q ~ b 2 q 
- M 2 )m q 4 


tf> = 


M a 2 q + b 2 q 
(mj - M 2 ) 2 2 ' 






g?< - 


= . 







(19) 

Here we take the zero quark mass limit. 

At loop level, the interactions of the WIMP with gluon are generated. The leading 
contribution is from one-loop diagrams including each heavy or light quark in Fig. |2J For 
the calculation of these diagrams, we adopt the Fock-Schwinger gauge for the gluon field: 

afAfa) = 0. (20) 

By using this gauge, we can express the gauge field by the field strength tensor as 

Afa) = \x P G a p »(0) + ■ ■ ■ , (21) 

so that we can extract out relevant effective interaction of the gluon easily. On the other 
hand, since the translation invariance is lost in this gauge, we have to distinguish two 
propagators for quarks in a gluon background, 



iS(p) = J d*xe^(T{q(x)q(0)})) , 

iS(p) = J d 4 xe-^(T{q(0)q(x)}) , (22) 



and those for squarks, 

zA(p) = f d 4 xe ipx (T{q(x)q j (0)}) , 

zA(p) = f d A xe- ipx (T{q(0)q\x)}) . (23) 

In the Appendix, propagators of the colored fermion and boson in the Fock-Schwinger 
gauge are collected for convenience. 

The two-point function of the neutralino in the gluon background (donated as T^) is 
calculated by using the above propagators; then the effective scalar coupling of neutralino 
and gluon is derived from the two-point function as 

f i -E*Mi (oa) 

JG\q — 5 l GG ' \ > 

where 

j^- [(a, - b ql5 )S(p + q)(a q + b ql5 )} A(g) 

7^4 [(o, + & g 7s)^(p - g)(a, - & g75 )] A(g) . (25) 

Here, the factor 1/2 comes from the neutralino being Majorana and "|gg" means to 
extract the coefficient for terms proportional to G a G a ^ u . For simplicity, we write T(p) to 
which the single flavor of quark (and squark) contributes. As we described in the previous 
section, we decompose it into short- and long-distance contributions and write as 

m , a„ fal + b 2 a 2 - b 2 n 

fb% = ^{^Mf + + <^m,f_y (26) 

The short- distance contribution comes from diagram (b) in Fig. [2} while the long-distance 
one is from diagram (d). The diagrams (a) and (c) vanish in our gauge. (See Appendix.) 
The mass functions are given as 



f- = m\B^ 

/i = ml(B^ + B[ 4 

t = B^+m 2 B^, (27) 



where 



d 4 q 



in 2 ((p + q) 2 - m 2 ) n (q 2 - rri^ m " 



g{n,m) 



ql \i q) 

74, 



^ ^ - V^ m) . (28) 



m 2 ((p + q) 2 — m^) n (q 2 — m 2 ~) m 




(a) 






(d) 



Figure 2: One- loop diagrams to generate interaction of the WIMP and gluon. Here, q 
and Q are light and heavy quarks, respectively, and q is for colored scalars, such as the 
squark in the MSSM. 



In the above expression, it is apparent that the typical momentum scale of integral is 
dominated by the squark mass in the short-distance contribution while it is dominated 
by quark mass in the long-distance one. With straightforward calculation of the integral, 
we have finally obtained 
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/i 
f- 

Here, 



(A — 6m~m 2 )(m~ + m 2 — M 2 ) 2m 2 mt 



g g 



6A 2 m? 



A 2 



3Am~ — (A — 6m?m 2 )(m? — m 2 + M 



A + \2m\m\ 
6A 2 



+ 



6A 2 m~ 

m^m 2 (m^ + m 2 — M 2 



+ 



rn 2 -rn 2 {m 2 ^ — m 2 — M 2 ) 

X 2 L ' 



A 2 



-L, 



3 Am 2 + 2(A + 3m 2 -m 2 )(m? 



m„ 



M 2 ) m 2 -(A + m 2 (m 2 --m 2 + M 2 )) 



6A 2 m 2 



M 4 - 2M 2 (m 2 q + m§) + (m 2 



A 2 



L . 

(29) 



m„ 



,2\2 



|A| 



. m 2 -+m 2 -M 2 + 

log — — 



i?+m,2_M 2 - 



|A| 



(A>0) 



2 tan" 1 



'|A| 



|A| 



m 2 +m 2 -M 2 
q q 



(A < 0) . 



(30) 



When we take zero quark mass limit, the mass functions are approximated as 



r~^ 



f S _ 

J+ 6m?(m 2 ~-M 2 ) 

2m? - M 2 



I s - 



6m? ( m ? - M 2 ) 2 
1 



6(m? - M 2 ) 2 ' 
j - 3m 2 (m?-M 2 ) ' l ' 



For the heavy quark loop in the long-distance contribution, it is checked that Eq. ([Tj 
is satisfied from Eqs. ( ITS]) and (115]) . just replacing q as Q, when m q — M ^> mg, as 
expected. The mass function f\ in the long-distance contribution is not singular in a 
limit of zero quark mass, while ft is. This is because f l + is proportional to m 2 before the 
loop momentum integral and the loop momentum integral leads to 1/m 2 

With all results, we obtain the effective coupling of the neutralino with gluon. As we 
described in the previous section, the contribution of light quarks should not be included in 
fa since the loop momentum has infrared cutoff, ~ GeV, implicitly. In Ref . [5] , however, 
the long-distance contribution of light quarks, i.e. fa\ q (q = u,d,s), is added in fa- 
This point can be seen by comparing our result with theirs in a simple expression under 
a certain limit. In the case of rriq ^> M, for instance, our result gives 

f ^ a s ^ ( a 2 q + b 2 q M\ a s y^ / a 2 Q + b 2 Q M a 2 Q - b 2 Q 1 \ 
J G ~^1^\ 24 ^1 J + 4^ ^ 1 24^^| + 12 m m 2 J ' ( ' 

g=all V '? / Q=c,b,t V 1 W q/ 

Focusing on the term which is proportional to a 2 + b 2 for the comparison, the result given 
in Ref. [5] has the form in which the term, —d s /4:Ti^2 uds (a 2 + 6 2 )M/24m|, is added 
to the above expression. We found that the long-distance contribution of light quarks, 
proportional to a 2 + b 2 , is regarded as the short- distance one in Ref. [5] and added in fa- 
Now we numerically show the cross section for the SI scattering of the bino-like neu- 
tralino with nucleon. The interaction of the bino-like neutralino with quark and squark 
is given as 

£ = y/2gY^2q(c^P R + c%P L )xQi + ^c- , (33) 

q,i 

where Pl/r = (1 =F 7s)/2 and gy is the gauge coupling constant of U(l)y. q± and q\ are 
the lighter and heavier squark, respectively, and the coefficients, c qL and c qR (i = 1,2), 
are 

C ql = - Y 1L COS H » C qR = Y QR Sm 9 q » 

cf L = Y qL sin 6 q , cf R = Y qR cos 6 q , (34) 
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Figure 3: x~P SI scattering cross section as a function of neutralino mass. Here we take 
parameters as (rriqL/M, m,q R /M, m q LR/M) = (1.5, 1.2, 0.1), and (2.0, 1.5, 1.0), and 
(3.0, 2.0, 1.0) from top to bottom. Solid lines show our result, while dashed-dot lines are 
the result in the case where light quark contribution in long-distance diagrams is included 
for reference. 



where Y UL = Y dh = 1/6, Y UR = 2/3 and Y dR = —1/3. The mixing angle of squark mass 

2m q fn qLR 



matrix 9 q is 



tan 29 f, 



m \h ~ m \ R ' 



(35) 



where we parametrize mass terms of left-handed and right-handed squark (q^ and g^) as 



a 



(mass) 



-m 



\lTlQl ~ m2 qR Q R QR + (m q rh qLR q* L q R + h.c.) . 



(36) 



From the interaction Lagrangian, a q and b q can be written as 

a q l) = —(-Y qL cos 9q + Y qR sin 9 q ) , b ( q l) = —(-Y qL cos 9 q - Y qR sin 9 q ) 



f/,. = —= (Y qL sin 9q + Y qR cos 9 q ) , b q 2) = — = (Y qL sin 9q - Y qR cos 9 q ) . (37) 



i (2) 
9 ^ 



V2" 



Then, we numerically calculate the SI cross section, using Eqs. ( [17]) and (TH?|) . (Here we 
include the QCD correction to the long-distance contribution to fa-) 

The cross section for SI scattering with proton, 0g D is plotted in Fig. [3] as a function 
of the neutralino mass. Here we give several results in the figure by taking (m q L/M, 
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rrig R /M, m qLR /M) = (1.5, 1.2, 0.1), (2.0, 1.5, 1.0), and (3.0, 2.0, 1.0) from top to bottom 
in solid lines. It is found that the cross section is suppressed as the neutralino mass 
become larger, and so it is unless masses of the neutralino and squarks are degenerate. 
In Fig. HI we present tree- and loop-level contributions to SI coupling in solid and dashed 
lines, respectively. The tree- level contribution, which mainly comes from g q in Eq. (fl9|) . 
is dominant, and the quark-squark loop contribution is subdominant; it accounts for 
O(10) % in the SI cross section. Just for reference, we also give the contribution from 
loop diagram in the case where the light quark contribution of long-distance is added in 
fa (dashed-dot lines in Figs. [3] and 0J. We have seen that the difference between our 
evaluation and the reference one is O(10) %. As shown in Fig. [3], the result is almost 
independent of the choice of the squark masses if they are in the same order of bino mass. 
This is because both fc D \q an d f(f\q & Ye the same order of contribution to fc in such 
a casejj The result, which has improved theoretical calculation of the SI cross section of 
dark matter with nucleon, will be important when we determine the local density and/or 
cross section of dark matter in the future experiments. Finally we also note that, in the 
public codes for studies of the neutralino dark matter, DarkSUSY [18] and MicroOMEGA 
[TU] . the gluon contribution to the spin- independent cross section is included by following 
Ref. [5]. Those programs should be corrected when one needs to calculate the cross section 
at less than O(10) % error. 

4 Wino-like neutralino dark matter 

In the anomaly mediation [20] the wino-like neutralino is a candidate for dark matter in 
the universe. The thermal relic abundance of the wino-like neutralino in the Universe is 
consistent with the WMAP observation when the wino mass is from 2.7 TeV to 3.0 TeV 
[2Tj . The lighter wino-like neutralino may be consistent with the observation of the dark 
matter abundance in the Universe if decay of gravitino or other quasistable particles may 
produce the dark matter nonthermally [22| [23]. 

The tree-level contribution to the cross section for elastic scattering of the wino-like 
neutralino with nucleon is evaluated in Ref. [21]. However, in the case that the superpar- 
ticles except for gauginos and the heavier Higgs bosons are much heavier than the weak 
scale in the MSSM, the tree-level interactions of the neutralino to quarks are quite sup- 
pressed. The scattering process is dominated by the W boson loop diagrams. Recently, 
we reevaluated the cross section for SI scattering of the wino-like neutralino with nucleon 
[8]. In the work, the two- loop contribution to the effective interaction of the neutralino 
with gluon is included, in addition to one-loop ones to the interaction with quarks. Then 
it is confirmed that both one- and two-loop contributions act as the leading contribution 
in the SI scattering process. Although calculation of two-loop diagrams is much involved, 



4 We note that there exist benchmark points to explain the observed dark matter abundance in which 
squark masses are much larger than the bino one. In such a case, the lightest neutralino should be a mixed 
state of bino and Higgsinos so that the SI cross section is determined by the Higgs boson exchange. We 
have checked that the long-distance contribution induced by the Higgs exchange is correctly evaluated in 
the previous works (e.g. Ref. [5]) and the numerical result of the cross section in such a case is unchanged. 
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Figure 4: Tree- and loop-level contributions in the effective coupling, f p . Here we take 
{rnqL/m^ mqR/m^, m q Lfi/m^) = (3.0, 2.0, 1.0). In the figure, tree- and loop-level ones 
are given in solid and dashed lines, respectively. Here, we also give a dashed-dot line 
for the case where contribution of light quarks in long-distance loops is added in fa for 
reference. 



it can be dealt with systematically in the Fock-Schwinger gauge. In the following, we 
show the technical detail of the calculation. 

Winos are gauginos which are partners of the SU(2)l gauge bosons in the MSSM, and 
they have the weak interaction. The relevant operators in our discussion is 



£ = -g 2 (xi fi x'Wl + h.c.) 



(38) 



where g% is the gauge coupling constant of SU(2)l- The wino-like neutralino accompanies 
the wino-like chargino (x~). The mass difference is dominated by one- loop contribution 
unless the Higgsino and wino masses are almost degenerate. 

Before getting down to the two-loop calculation, we quickly review one-loop contri- 
bution to the SI interaction of x-iV '. The effective interactions of x with light quarks in 
the effective Lagrangian are generated by diagrams in Fig. |5j The diagram (a), which is 
induced by the standard-model Higgs boson (h°) exchange, contributes to f q , while the 
diagram (b) generates the other terms in the effective Lagrangian. The results are given 
in Ref. [8]. The obtained effective couplings are finite in the limit M — > oo. Thus, the SI 
interaction of x-iV is not suppressed even if the wino-like neutralino is much larger than 
the W boson mass, as pointed out in Ref. [TJ. 

Next, let us derive the effective interaction of the neutralino with gluon, which is 
generated at two-loop level. Three types of diagrams in Fig. [6] contribute to fa- Diagram 
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Figure 5: One-loop contributions to effective interactions of the wino-like neutralino and 
light quarks. 




(a) 





(c) 



Figure 6: Two-loop contributions to interactions of wino-like neutralino and gluon. 



(a) includes heavy quark loops, and it is evaluated from fq as given in Eq. ffTBT) . On 
the other hand, we need to calculate irreducible two-loop diagrams (b) and (c) explicitly. 
In order to handle this calculation, it is convenient to evaluate the vacuum polarization 
tensor of W boson in the gluon background, 



*lWg) 



-E 

CM) 



d 4 p gj 

(2tt) 4 2 



Tr 



li,PLSd{p)PRluS u {P - q) 



(39) 



Here, S u and Sd are propagators of up-type and down-type quarks, respectively, in the 
Fock-Schwinger gauge for gluon field. The vacuum polarization tensor is decomposed to 
the transverse and longitudinal parts as 



IV(g) = U T (q 2 )(-g^ + 



Wv 



+ n L ( g : 



2x %<lv 



(40) 



q" q* 

We found from an explicit calculation that the longitudinal part of the self energy does 
not contribute to fc- Thus, what we have to calculate is coefficients for terms with gluon 
scalar operator G a tlv G aiiV in Il T (g 2 ), i.e. Il T (g 2 )| GG . 
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For the calculation Il r (g 2 )| GG , we furthermore decompose the contributions from di- 
agram (b) and (c) in Fig. [H]as follows; 



n£V;0loo 



a 2 a s 

~2A 

a 2 a s 

~T 

a 2 a s 



6B. 



(2,2) 



,(2,2) 



20 



+ Q -°22 + Q B \ 



(2,2) 



-m„. 



-m„ 



-2B. 



-2B 



(4,1) 
20 

(1,4) 
20 



q 2 B : 



22 



fB^ 






(41) 



where i is the flavor index and a 2 = g 2 /A-K. Tv T ' is obtained from diagram (b) and 11^ 
(Ilj? ) is from diagram (c) in which gluon lines are attached with the propagator of up- 
(down-) type quark. In the above expressions, loop integrals are defined as 



d A q 



% 



in 2 ((p + q) 2 — m 2 1 ) n {q 2 — m^) r 
dS 



%<lv 



in 2 ((p + q) 2 — m 2 ) n (q 2 — m 2 



T->(n,m) 



-°20 Utiv T J-> 22 PfiPu 



(42) 



where m u and m^ are masses of up-type and down-type quarks, respectively. 

Notice that, in the loop integrals for the self-energy part of diagram (c), the typical 
momentum is the mass of quark which emit two gluons, i.e. the up- and down-type quark 
mass in W£ and ITj? , respectively. This fact means that diagram (c) gives the long- 
distance contribution in the part of one-loop self-energy of the W boson. Therefore, we 
add the contribution of lij, for i = 2 and 3 (up-type quark is heavy one) and 11^ for 
i = 3 (down- type quark is heavy one). On the other hand, the loop momentum of quark 
loop in diagram (b) is dominated by the external momentum of the quark loop diagram; 
therefore, all quarks contribute in the loop. Then Ht{q 2 )\gg is properly obtained as 



iW)| GG = £ 4V;*)| G c+£4 cl V 



| GG + np> 



(s 2 ;3)| 



GG 



(43) 



=1,2,3 



=2,3 



Executing the integral and neglecting quark masses except for the top one, we finally get 



n T ( g 2 )| 



GG 



ct s a 2 



1 2q 2 



3m 2 



24 (q 2 



m 



2\2 



1 1 

+ 2 x 

12 q 2 



(44) 



The first term in the right-handed side comes from loops of the third generation, while 
the second term is induced by the first- and second-generation quarks. We observed that 
the singular terms in 1/m 2 are cancelled in the loop integrals of 11^ ' c so that 11^ ' c 
vanishes in the m q — > limit. Actually, this is expected from the indication of Eq. (j!4p 
and the fact that diagram (b) in Fig. [5] does not contribute to the neutralino-quark scalar 
coupling f q in a limit of zero quark mass [8] . 
Combining the above results, we get fa as 
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Figure 7: \~V SI scattering cross section as a function of neutralino mass M for wino-like 
neutralino case. Here we take the Higgs boson mass as 115, 130, 300 GeV, and 1 TeV from 
top to bottom. Results where the QCD correction is taken into account for long-distance 
contribution are given in solid lines, and those without the QCD correction are in dashed 
lines. 



where m^o and mw are the Higgs and H^-boson masses, respectively, and y = mf/M 2 . 
In addition, we include the QCD correction to the long-distance contribution of heavy 
quarks to fa- The mass functions ^(x), gB3(x,y) and ^bi(x) are given in Ref. [7|. 

In Fig. [7J we show the numerical result for the SI cross section. Here the evaluation 
with (without) the QCD correction is given by taking m h o = 115, 130, 300 GeV, and 
1 TeV. While the latter two values may not be realistic in the MSSM, the next-minimal 
supersymmetric standard model, for example, may predict a larger Higgs boson mass. It 
is found from the figure that the cross section is more suppressed unless the Higgs boson 
mass is larger than a few hundred GeV. This fact is the consequence of an accidental 
cancellation in SI coupling as it is shown in Ref. [8]; the QCD correction makes the 
contribution from the Higgs boson exchange larger so that the cancellation works hard|j 
On the other hand, the effect of the cancellation becomes weaker as the Higgs boson mass 
is larger. 



5 In Ref. [8] , the QCD correction for long-distance contribution and the contribution from the twist- 2 
operators of bottom and charm quarks are not included. 
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5 Conclusion 

In this article, we have systematically calculated the spin-independent cross section of the 
nucleon-lightest neutralino scattering process at loop level. Although the one-loop (higher- 
loop) diagrams are usually the higher-order of quantum correction to tree-level (one- 
loop) ones, such higher order diagrams of neutralino-gluon scattering work as the leading 
contribution. We have evaluated the effective coupling of neutralino-gluon scattering by 
explicit calculation of loop diagrams in the Fock-Schwinger gauge. In this gauge, we can 
calculate systematically the loop diagrams, tracking the loop momentum scales which 
dominate in the integrals. Making the most of it, we have separately evaluated short- and 
long-distance contributions in the loop diagrams, i.e. the contributions in which the loop 
momentum is characterized by masses of new particles heavier than quarks (like mass 
scale of neutralino or squark) and those of quarks. Being aware that the contribution 
of light quarks must not be included in the long-distance diagrams, we have numerically 
calculated the spin-independent effective coupling to give the cross section of nucleon- 
neutralino scattering. In order to focus on the contribution of neutralino-gluon scattering, 
we considered the bino- and wino-like neutralino cases. In the former case, we have found 
that the effective coupling of neutralino-gluon scattering differs by 20-30 % from the one 
evaluated in the case where light quarks are included in the long-distance contribution 
incorrectly; it gives O(10) % alteration in the SI cross section. This improvement in the 
theoretical calculation of the cross section will be essential when we explore the nature (i.e. 
interaction and/or local density) of dark matter in the future experiments. On the other 
hand, in the latter case, the quark- and gluon-neutralino scattering contributes in almost 
the same magnitude and opposite sign. Then, they cancel each other to suppress the 
SI cross section. Here we have also taken into account QCD correction for long-distance 
contribution; consequently, we have found that the cancellation works hard. 

Acknowledgment 

The work was supported in part by the Grant-in-Aid for the Ministry of Education, 
Culture, Sports, Science, and Technology, Government of Japan, No. 20244037, No. 
2054252 and No. 2244021 (J.H.) and Research Fellowships of the Japan Society for the 
Promotion of Science for Young Scientists (K.I.). The work of J.H. is also supported 
by the World Premier International Research Center Initiative (WPI Initiative), MEXT, 
Japan. 

A Propagators in Fock-Schwinger gauge 

Here we present colored fermion and scalar boson propagators in the Fock-Schwinger 
gauge for gluon, and useful formulas to evaluate effective coupling of WIMP and gluon. 
We also show derivation of the coupling in a model in which the WIMP has an interaction 
with the Higgs boson, as a simple exercise. 
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Figure 8: Fermion propagators. 



Defining gluon field A^ in covariant derivative, D^ = <9 M + ig s A a ^T a , the gauge- fixing 
condition in the Fock-Schwinger gauge is written as 



x»A°(x) = . 



(46) 



In this gauge, the vector field can be directly expressed in the terms of the gluon field 
strength tensor as 



Afa) = da ax p G a pfl (ax) 

Jo 

= ^^(0) + ^,V(D a G w (0)) a + ^V/(D a i}^(0))" + ■ ■ • . 

(47) 

While the gauge-fixing condition breaks the translation invariance, it is recovered in the 
gauge-invariant objects. The propagators of colored fermion in the gluon background are 
given as 
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iS(p) = / d A x e ipx (T{tP(x)iP(0)}) 
= tS {0 \p) 



+ Jd% iS®(p) g sl a Q G ^^ (4) (^i)) iS^{p- h) 
+ Jd%d% iS®(p) g s i a Qg^ JUwfo)) iS(0) (P - fc i) 9s^ 
x (l G e^ 8i4) ( k *)) iS(0) (p ~h-k 2 ) + --- , (48) 

iSip) = f d 4 xe- ipx (T{tP(0)tP(x)}) 
= tS {0) (p) 
+ Jd% iS^ip + h) g sl a (hs a ^8^{p)\ iS^ip) 



+ Jd%d 4 k 2 iS<®(p + h + k 2 ) g sl a {\ G ^^^\k 2 ] 

xiS®(p + h) gtf flgG^-^-S^ih^j iS®(p) + • • • , (49) 

where iS^°'(p) = i/(p ' — m) and G^ v = G a ^JT a . In the actual calculation, terms including 
covariant derivatives, such as D a G^ y and D a D^G^ u , are ignored, since they are irrelevant 
to evaluation of the SI cross section. In text we evaluate the scalar coupling of WIMP 
with gluon, fa, by extracting out the bilinear term of the gluon field strength from the 
WIMP (W boson) two-point function in the bino- (wino-)like neutralino dark matter case. 
The scalar operator of gluon G p ] v G a ^ 1 ' is projected out from the bilinear term of the gluon 
field strength as 

GanGpv — -^p a P aG ap(7 {g^g^y - g av gpy) 

~ 2 9a ^°^ ~~ 2 9fiu ° 9al3 ~~ 2 9au ° 9 ^ ~ 2 9 ^° 9av 

+ O^v, (50) 

where 9 pv is the twist-2 operator of gluon (Eq. (Tjfl)) and O a ^ v is a higher- twist operator, 

sr\g __ /~ia /~ia 

1 1 1 1 

- 7^9apG p G pu — -g^vG a G p/3 + -g au G pG pp + -gp^G a G pv 

+ l -G a pa G a ^(g aP g, v -g au g p ,). (51) 
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We also give colored scalar propagators: 

»A (P ) , /^<r WW t ( o, } ) 

= iA i0 \p) 

+ J d 4 h iA<°>(p) g s (2p - fcO a (^^^"^(fci)) ^ A(0) (P - *i) 

+ Jd 4 hd 4 k 2 iA(°)(p) <? s (2p - feO a Q^^-^i)) ^A (0) (P " fei) 

x^ s (2p - 2fci - fc 2 )" Qg^^U< 4 >(A; 2 )) iA^(p - fci - fe 2 ) 

+ jd 4 k 1 d 4 k 2 iA^\p){-igl){^-G a ,^ 4 \k 1 )^ {^G\^-f\k 2 ) 

xiA i °\p-k 1 -k 2 ), (52) 



iA(p) = / d 4 :r e^T^O)^)}) 
= iA (0) (p) 



+ y^ iA^ip + h) 9s (2p + k 1 ) a 0^G a ,^-5^ A \h)\iA^\p) 
+ /" d%d% iA<°>(p + fe! + fc 2 ) ^ s (2j9 + 2A;i + A; 2 ) a Qg qm ^|-<5( 4 )(A; : 



2, 
9 



xiA^ip + h) g s (2p + k 1 f r-G^^-d^ikM iA^ip) 



I d 4 k 1 d 4 k 2 iA^°\p + k 1 + k 2 



x <-'^» O-d^^'V 0w <4)(fe) ) !Al0,(p) • (53) 

where iA^\p) = i/(p 2 — m 2 ). The scalar propagators are reduced in a more convenient 
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Figure 9: Scalar boson propagators. 



form for practical usage as 
iA(p) = iA {0 \p) 



/ ( /7., ( J-SMfa)) gp a G a ,{iA^{p)f 



dk 



in 



I f/V/'A-' ( J^ (4) (^i)l f^ (4) (^) ) sW^G^A^)) 3 



dA; 



Ip 



0fc. 



2u 



(fkn 



d 



dk 



2// 



5^(k 2 ))g 2 p a G a ,G^A^(p)f 



+ Jd%d%iA^{p)(-ig 2 )f^G 











,,jjj~5 {A \h))[^G a u ^-5^(k 2 , 



2 v dk 2v 



xiA ( -°\p-k 1 -k 2 ) 

iA(p) = iA (0 \p) 



(54) 



+ / d% 







6^(k 1 ))gp a G a JiA^{p)f 



dk lfl 
I <l l ^<l% (J-SW(kS) (4~^\k2) ) g 2 p a P P G a ,G^A^(p)) 3 



dk 



ip 



0fc. 



2i> 



9 



/ d% ( ^<$ (4) (A; 2 ) gYG^G^iA^ip)) 3 



dk 



2i/ 
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(55) 
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(b) 



Figure 10: Effective interactions of \ with light quarks and gluon, which are induced by 
Higgs boson exchange. 



In the calculation of the SI cross section of the bino-like neutralino in Sec. [31 the diagrams 
(a) and (c) in Fig. [2] are found to be zero in the Fock-Schwinger gauge, since terms except 
the first and fifth terms in iA(p) and iA(p) vanish in the calculation. 

Last, we show the loop diagram calculation in the Fock-Schwinger gauge for a simple 
scenario as an exercise. We consider a model in which the WIMP x interacts with the 
standard-model Higgs boson; 



C 



r,c,uO 



-c XXh 



(56) 



Effective interactions of x with light (and heavy) quarks generated by diagram (a) in 

Fig. ma 



fq(Q) 



cg 2 



2m w m 2 h0 



(57) 



at leading order. Other coefficients in Eqs. (EJ) and (EI) are zero at the leading order. The 
effective coupling of WIMP with gluon is induced by diagram (b) in Fig. HOj 



./. 



G 



- E f^ j wf l[lS{p)] 



GG 



Q=c,b,t 



Air 



v^ , f d A p m Q p 2 

h: QmQ J ^(p 2 -m 2 Q r 



Q=c,b,t 

12vr X J<3 ' 



(58) 



Here we neglected the QCD correction just for simplicity. 
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